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Within a dynamical quark recombination model we explore various proposed event-by-event ob- 
servables sensitive to the microscopic structure of the QCD-matter created at RHIC energies. Charge 
fluctuations, charge transfer fluctuations and baryon-strangeness correlations are computed from a 
sample of central Au-|-Au events at the highest RHIC energy available (-ysjvjv=200 GeV). We find 
that for all explored observables, the calculations yield the values predicted for a quark-gluon plasma 
only at early times of the evolution, whereas the final state approaches the values expected for a 
hadronic gas. We argue that the recombination-like hadronization process itself is responsible for 
the disappearance of the predicted deconfinement signatures. This might explain why no fluctu- 
ation signatures for the transition between quark and hadronic matter was ever observed in the 
experimental data up to now. However, it might also be interpreted as a clear indication for a 
recombination like hadronization process at RHIC. 
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It is widely believed that a phase transition from a 
quark-gluon plasma (QGP) to hadronic matter occurs in 
central ultra-relativistic heavy- ions coUisions at RHIC. 
In order to study the properties of the extremely heated 
and compressed matter created in these events, numer- 
Drobes based on fluctuations have been proposed 
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iBTlllIjT For a compre 



hensive overview in the physics of event-by-event fluc- 
tuations we refer the reader to [2^. Among them, es- 
pecially charge ratio fluctuations, charge transfer fluctu- 
ations and baryon-strangeness correlations were promi- 
nently proposed to pin down the formation of a decon- 
fined phase at RHIC ij, IM HI, IH [H, 113, [Sll, IH, [H . 

These observables are based on event- by- event fluctu- 
ations of conserved charges within a given rapidity range 
and are sensitive to the microscopic nature of the matter. 
It was pointed out that these quantities reflect the prop- 
erties of the system in the first instant of the collision 
and should survive the whole course of the evolution of 
the system. The argument in favour of the survival of the 
signal through all the stages of the collision for the above 
mentioned fluctuations probes is the following: With a 
strong transverse and longitudinal flow, locally conserved 
quantities (charge, baryon number and strangeness) will 
be frozen in a given rapidity window because the expan- 
sion is too quick for the charges to move out of the respec- 
tive rapidity slice. Thus, if a QGP is created, the fluctu- 
ation of these quantities should survive further evolution 
through the hadronic phase. 
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It is clear that the size of the rapidity window for the 
fluctuation study must not be too wide in order to avoid 
global conservation which would lead to a vanishing sig- 
nal, but also neither too small to avoid purely statistical 
fluctuations and the transport of charges in and out of 
the window by hadronic rescattering. The generally ac- 
cepted rapidity width is of the order of Ay = 0.5— 1 units 
in rapidity. In contrast to the RHIC energies explored 
here, it was argued that the diffusion rate for secondaries 
at the CERN-SPS might be strong enough to blur the 
fluctuation signal almost to the (observed) resonance gas 
value [H. 

A key point that is usually not addressed in the discus- 
sion of fluctuation signals is the influence of hadroniza- 
tion itself. A possible mechanism for the parton-hadron 
transition is the recombination of quarks and anti-quarks 
into hadrons [H, [H, [13, [H, HE] ■ Elliptic flow and nuclear 
modification factor Raa measurements at RHIC [10, |4l| 
have given strong evidence supporting recombination as 
the mechanism responsible for hadronization. A first ex- 
ploratory study on the influence of parton recombination 
on charge fluctuation was performed in |4^ . There it 
was shown that the coalescence of quarks through the re- 
combination mechanism does indeed lead to results com- 
patible with the available experimental data on charge 
fluctuations at RHIC. 

In this paper, we study charge ratio fluctuations, 
charge transfer fluctuations and baryon-strangeness cor- 
relations with a dynamical recombination model (the 
quark Molecular Dynamics model, qMD [H, 0, H^] ). 
To pin down the influence of the hadronization process in 
detail we explore the suggested quantities over the whole 
time evolution of the system from the pure quark stage 
to the final hadrons. The set of events consist of central 
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Au+Au collisions at the highest RHIC energy available 
{^/sNN = 200 GeV). We will finally conclude that the 
hadronization process itself is responsible for the change 
of all investigated observablcs from the initially partonic 
value to the finally observable hadronic value. Thus, pro- 
viding evidence for a recombination like hadronization 
mechanism at RHIC energies. 

The qMD model [H, H [H 13 employed here 
is a semi-classical molecular dynamics approach where 
quarks are treated as point-like particles carrying color 
charges and interact via a linear heavy quark potential. 
Initial conditions [ss] for the qMD are taken from the 
hadron-string transport model UrQMD [13, After 
the two incoming nuclei have passed through each other, 
(pre-)hadrons from the string and hadron dynamics of 
the UrQMD model are decomposed into quarks with cur- 
rent masses m„ — nid = 10 MeV and ~ 150 MeV. 
At the highest RHIC energy, this happens at a center of 
mass time oi t = 0.15 fm/c. The quarks are then let to 
evolve and interact within the qMD via a linear poten- 
tial F ( | — | ) = K I Ti — r j I , where k is the string tension 
and r„ is the position of particle n. Therefore the full 
Hamiltonian of the model reads: 

" / 1 

1=1 ij 

Where N counts the number of particles in the system 
and the term takes into account the color dependence 
of the interaction. 

The quark-(anti-)quark interaction within this poten- 
tial naturally leads to confinement through the binding 
of (anti-)quarks into color neutral clusters. New hadrons 
are formed from quarks whose momentum and position 
are close to each others. Typical values for the relative 
momenta of the quarks in the two-particle rest frame at 
hadronization are \pq\ — \pq\ < 500 MeV, the typical dis- 
tance is below 1 fm. Hadronization thus occurs locally 
into hadronic clusters of mesonic and baryonic type that 
resemble the Yo-Yo states of the LUND model. These 
clusters are allowed to decay in the further evolution of 
the system and the hadronization process therefore al- 
lows to conserve entropy. 

As the initial state of the system is color neutral, all 
the quarks of the system will eventually gather into color 
neutral clusters. Electric charge and strangeness are con- 
served during the whole evolution of the system, i.e. 
s — s = and the net-charge equals the initial ch arg e of 
the incoming nuclei. The reader is referred to (isl. l44 [isjl 
for a detailed discussion of the qMD model. Note that in 
the present calculations, only u, d and s quarks are in- 
cluded. Furthermore, all parton production occurs in the 
early stage of the reaction during the UrQMD evolution. 
There is presently no mechanism to create new (di-) quark 
pairs during the qMD evolution stage. Thus, the present 
model provides an explicit recombination transition from 
quark matter to hadronic matter in a dynamical and ex- 
panding medium. 




time [fm] 

FIG. 1; Fraction of the number of quarks from the total 
number of particles as a function of time at midrapidity for 
Au+Au reactions at y^sjvjv = 200 GeV. The arrow depicts 
the hadronization time. 

Let us set the stage by exploring the time evolution 
of the hadronization dynamics in the model. Fig. [T] de- 
picts the fraction of quark matter on the total num- 
ber of particles in the system (i.e. quark fraction = 
(uq + n-q) I (nhadron + "-g + "-q)) as a functiou of time. One 
observes that the fireball stays in a deconfined state dur- 
ing the first 6 fm/c where almost no quarks hadronize. 
As the system expands further and the density decreases, 
quark recombination into baryons and mesons occurs and 
the number of deconfined quarks drops to zero. Next 
we turn to the investigation of the various fluctuation 
signals. The electric charge ratio fluctuations were pro- 
posed as a clear signal for the onset of the quark-gluon 
plasma phase [29^1 . The basis for the argument is that the 
quanta of the electric charge are smaller in a quark gluon 
plasma phase than in a hadron gas and are distributed 
over a larger number of particles. Moving one charged 
particle from/to the rapidity window then leads to larger 
fluctuations in a hadron gas than in a QGP. The electric 
charge ratio fluctuation can be quantified by the measure 
D defined as: 

D = ^{N,n)m^)^y ■ (2) 

Where N^h stands for the number of charged particles, 
R = (l-hF)/(l-F) withi^ = Q/Vch, Q being the electric 
charge. Following [l^, charge fluctuations are corrected 
with the factors and Cy to take into account the fi- 
nite acceptance. As suggested in [H,!!^, the quantity D 
is calculated in a rapidity window of y = ±0.5. It was 
argued that depending on the initial nature of the sys- 
tem, D will yield distinctly different results: = 1 for 
a quark-gluon plasma, D = 2.8 for a resonance gas and 
_D = 4 for an uncorrelated pion gas. 
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FIG. 2: Corrected charge fluctuations _D as a function of 
time within the qMD model for Au+Au reaction at ^sjvjv ~ 
200 GeV (full symbols). Also shown are the values for an 
uncorrelated pion gas, a resonance gas and a quark-gluon 
plasma. The arrow depicts the hadronization time. 



FIG. 3: Charge transfer fluctuations at midrapidity and — 
as a function of time within the qMD model (full symbols) for 
Au+Au reactions at ^sjvjv ~ 200 GeV. The arrow indicates 
the hadronization time. 



Experimentally, charge ratio fluctuations have been 
measured at RHIC energies by STAR ^ [s^l and 
PHENIX [51I, [l^l . Both experimental analyses yield re- 
sult compatible with a hadron gas - in strong contrast to 
the first expectations. Further results from the CERN- 
SPS mill based on a slightly different measure for the 
charge ratio fluctuations did also yield results compat- 
ible with the hadronic expectation. Fig. [2] shows the 
result for D from the qMD recombination approach as a 
function of time. In the early stage, when the system is 
completely in the deconfined phase, Z? = 1 as expected. 
When approaching the hadronization time, D starts to 
increase and reaches I? « 3.5 after hadronization. As 
can be seen from Fig. [T] the increase of D occurs exactly 
at the same time as the recombination of the quarks and 
anti-quarks to hadrons proceeds. The slight decrease of 
D at later times is related to the decay of resonances. 

As a next observable, we now turn to charge transfer 
fluctuations that were also suggested to provide insight 
about the formation of a QGP phase. Charge transfer 
fluctuations are a measure of the local charge correlation 
length. They are defined as [III, HI]: 



(3) 



with the charge transfer uirj) being the forward-backward 
charge difference: 



"('7) = [Qpiv) 



(4) 



where Qp and Qb arc the charges in the forward and 
backward hemisphere of the region separated at 77 = 0. 
In our calculations, we take a total window of y = ±1, 



corresponding to the STAR acceptance. Experimental 
data on this oberservable is not available up to now. 

Because the measured quantity is local, it can give in- 
formation about the presence and the extent of a QGP in 
rapidity space. Thus, one expects to observe the lowest 
value of the charge transfer fluctuations at midrapidity, 
where the energy density is the highest and where the 
plasma is located. The local charge fluctuation is ex- 
pected to be much lower in a quark-gluon plasma than 
in a hadron gas. The results from the present calculations 
are shown in Fig. [3] As expected, the correlation length 
(at central rapidities) is small, with Du/ (dNch/dy) « 0.1, 
as long as the system is in the quark phase. However, 
similar to the charge ratio fluctuation discussed above, 
the charge transfer measure increases with time up to its 
hadronic value of Du/{dNch/dy) « 0.5 when hadroniza- 
tion has happened. The final state result is in agreement 
with the value given by HIJING calculations and there- 
fore in line with the hadronic expectation [5^ . 

Finally, we analyse the baryon-strangeness correlation 
Cbs 1231 ■ This correlation was proposed as a tool to 
study the property of the matter created in heavy ion 
collisions. The baryon-strangeness correlation is defined 
as: 



Cbs — —3 



{BS)-{B){S) 
(52) - (5)2 



(5) 



where B and S are the baryon number and strangeness 
in a given event [soj . 

The rationale behind this quantity is the fact that 
baryon number and strangeness are differently correlated, 
depending on the phase the system is in. In an ideal 
weakly coupled quark-gluon plasma, strangeness will be 
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FIG. 4: Cbs correlation coefficient as a function of time for 
central Au+Au reactions at ^/snn = 200 GeV (full circles). 
Also shown are the predicted values for a hadron gas and a 
QGP. The arrow depicts the hadronization time. 

carried by strange quarks and is therefore strictly cou- 
pled to baryon charge. Thus, a clear correlation between 
baryon charge and strangeness is expected in a quark- 
gluon plasma. The expected numerical value for an ideal 
QGP is Cbs = 1 [H, H^l- In a hadron gas on the con- 
trary, strangeness can be carried without baryon number 
(e.g. with strange mesons). As a result the correlation 
between strangeness and baryon number will be weak- 
ened compared to the quark matter scenario. The numer- 
ical value for a non-interacting hadron gas is Cbs = 0.66 

MM- 

The behaviour of Cbs as a function of time for the 
dynamical recombination model under study is depicted 
in Fig.m For early times, Cbs starts from the expected 
value of unity in agreement with the ideal weakly coupled 
quark-gluon plasma value. During the course of the re- 
combination of the quarks, Cbs approaches the hadron- 
gas value Cbs ~ 0.6. 

In conclusion, we have studied a variety of suggested 
event-by-event signatures for the formation of a decon- 



fined QGP state within a dynamical quark recombination 
approach. 

The analyses was done for central Au+Au events at 
y'SNN = 200 GeV and involved charge ratio fluctuations, 
charge transfer fluctuations and baryon-strangeness cor- 
relations. For all these predicted "smoking gun" QGP 
observables, we find that the hadronization by recombi- 
nation leads to results expected for a hadron gas in the 
final state. This is especially remarkable, as the initial 
values for these observables were identical to the pre- 
dicted QGP values. 

For all these quantities, the change of the observables 
from there QGP value to the hadronic gas value can be 
traced back to the recombination hadronization mecha- 
nism because the change of the quantitative values of D, 
Du and Cbs takes place during the time of hadroniza- 
tion. 

From these observations we draw two mutually con- 
verse conclusions: 

1. The influence of the recombination/hadronization 
on fluctuation probes is strong enough to blur the 
initially present QGP signature. This might ex- 
plain why fluctuation measurements have not pro- 
vided the expected proof for the formation of a 
plasma of quarks and gluons. 

2. However, if one assumes that a QGP was indeed 
formed at RHIC energies, the experimental fact 
that all of the discussed fluctuation probes turn 
out to yield the hadronic value can be seen as a 
strong argument supporting recombination as the 
mechanism responsible for hadronization at RHIC. 
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It should be noted that the qualitative results of the 
present study are not restricted to any specific initial 
state. The UrQMD model is solely used to provide an 
exemplary initial state after the initial qq production has 
taken place. 

In Ref. [5^] it was shown that the true value of Cbs 
can be reasonably approximated by taking into account 
charged kaons and A's only. This might be easier for cur- 
rent day experiments because it avoids the measurement 
of neutral particles. 



